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Abstract 

Parton distributions in a valon in the next-to-leading order is used 
to determine the patron distributions in pion and kaon. The validity 
of the valon model is tested and shown that the partonic content of 
the valon is universal and independent of the valon type. We have 
evaluated the valon distribution in pion and kaon, and in particular it 
is shown that the results are in good agreement with the experimental 
data on pion structure in a wide range of x — [1CP 4 , 1]. 
PACSnumbers: 14.40.-n, 12.39.-x, 14.65.-q 



1 INTRODUCTION 

Unlike the structure function of proton, there are relatively fewer informa- 
tion on the structure of meson and in particular pion. The parton distribu- 
tions in proton have been studied extensively, both theoretically and experi- 
mentally, in a wide range of Q 2 = [0.45, 10000] GeV 2 and x = [10~ 5 , 1], Yet, 
pion plays an important role in QCD and its presence is felt everywhere in 
hadron physics: pion cloud of nucleon, baryon-meson fluctuation, and decay 
of quark to pion-quark which can explain certain aspects of flavor symmetry 
breaking in the nucleon sea, are just a few to name . As a result it is impor- 
tant to determine and understand its internal structure, which also renders 
useful information about nonperturbative QCD. However, at present, the 
parton distribution functions of pion are far from being satisfactory. 
The meson structure is measured in a number of Drell-Yan processes [1] [2]. 
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Such measurements are concentrated in the intermediate and large x region, 
mostly above 0.2, and hence, pertinent to the valence quark distribution. 
Recently, ZEUS Collaboration have published pion structure function data 
at very low x values from the leading neutron production in e + p collision 
[3] which provides some details about the sea quark distribution in pion. 
Unfortunately, there exists ambiguity in the normalization of ZEUS data. 
ZEUS collaboration have used two different methods to normalize the data. 
The results differ by a factor of two. In a recent paper [4] attempts are made 
to clarify the normalization ambiguity in the ZEUS data and independently 
calculate the pion structure function based on the valon model. The valon 
model essentially treats the hadron as the bound state of its valons. Each 
valon has its own partonic structure calculated in QCD. Measurements of 
Natchmann moments of proton structure function at Jefferson Laboratory 
[5] makes the valon model more credible. The findings of Ref.[5] points to 
a new type of scaling which can be interpreted as a constituent form fac- 
tor, consistent with the elastic nucleon data. This, in turn, suggests that 
the proton structure originates form the elastic coupling with the extended 
objects inside the proton [6]. If confirmed, such an extended object can be 
identified as valon. 

Of course, the notion of structureful objects in hadrons is not new. Altarelli 
and Cabibo [7] have used the concept in the context of SU(3) x 0(3) and 
R.C. Hwa has termed them as valons and further developed the concept 
and showed its application to many physical processes [8]. It is now well 
established that one can perturbatively dress a QCD Lagrangian field to all 
orders and construct a structureful object (valon) in conformity with the 
color confinement [9] [10]. More recently, the partonic content of a valon 
is calculated in the Next-to-Leading Order (NLO) [11] and shown that if 
convoluted with the valon distribution in proton, it gives a fairly accurate 
description of the proton structure function data in the entire kinematical 
range of measured values. 

The underlying assumption of the valon model is that the structure of a 
valon is independent of valon type and the hosting hadron. Therefore, it 
should provide insight into the structure of hadrons other than proton, for 
which the experimental information are either a rarity or less accurate. In 
fact, in Ref. [4] the method described in Ref. [11] is used to calculate the 
pion structure function. However, the focus was on the low x region and the 
knowledge gained was pertinent to the sea quark distribution in the pion; 
and to provide a resolution of the normalization ambiguity. Thus, the aim 
of this paper is three-fold: (a) to extend the analysis to include the inter- 
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mediate and large x region; (b) To test the validity of the valon model and 
extract the valon distribution of mesons, and (c) to elaborate on the kaon 
structure. 

2 The Valon Model 

In the valon model, we assume that baryons and mesons consist of three, 
and two valons, respectively. Each valon contains a valence quark of the 
same flavor as the valon itself and a sea of partons (quarks, antiquarks, 
and gluons). At low enough Q 2 the structure of a valon cannot be resolved 
and the hadron is viewed as the bound state of its valons. At high Q 2 the 
structure of a valon is described in terms of its partonic content. For a 
U-type valon, say, we may write its structure function as 

F?(z, Q 2 ) = \z[G± +Gu) + l -z{Ga + G| + G± + Gs_) + ... (1) 

where all the functions on the right-hand side are the probability functions 
for quarks having momentum fraction z of a U-type valon at Q 2 . A similar 
expression can be written for other types. In Ref. [11] the probability 
functions, or parton distributions in the valon are calculated in QCD to 
the NLO at the scales of Ql = 0.283 GeV 2 and A = 0.22 GeV. Without 
going into the details, it suffices here to give the functional forms of them 
as follows, 

zq v v %™ ce (z,Q 2 ) = az b (l-z) c , (2) 

zq^T(z, Q 2 ) = az\\ - z?\\ + V z + ^°- 5 ]. (3) 

The parameters a, 6, c, a, (3, 7, r], and £ are functions of Q 2 and are given in 
the appendix of Ref. [11] . Gluon distribution in a valon has an identical form 
as in Eq. (3) but with different parameter values. Eqs. (l)-(3) completely 
determine the partonic structure of the valon without any new parameter. 
We note that the sum rule reflecting the fact that each valon contains only 
one valence quark is satisfied for all Q 2 : 

[ 1 q V v Tnce(z,Q 2 )dz = l. (4) 
JO 
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3 The Meson Structure Functions 



3.1 A. Pion 

The determination of parton content of hadron requires the knowledge of the 
valon distribution in that hadron. Let us denote the probability of finding 
a valon carrying momentum fraction y of the hadron by Gvaion(y), which 

h 

describes the wave function of hadron in the valon representation, containing 
all the complications due to confinement. Following [4, 8, 11] and [12], we 
write the valon distribution in a meson as: 

G vaion (y) = — 1 -y m (\-y) Vm . (5) 

meson P[fJ-m + 1 5 + 1] 

with the requirements that the above form satisfies the number and momen- 
tum sum rules: 



I 
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G valon dlj = 1 / llG valon d,V = 1 (6) 

meson J Q meson 

valon 

where, (3[i,j] is the Euler beta function and G vaion (y) stands for the distri- 

h 

bution of a U- valon in ir + or a D-valon in ir~ . By interchange of \x <-> v the 
anti-valon distribution in the same meson is obtained. 

An essential property of the valon model is that the structure of hadron in 
the valon representation is independent of the probe. This means that the 
parton distribution in a hadron can be written as the convolution of the 
partons in the valon and the valon distribution in the hadron. For the case 
of pion, this translates into: 



XUvalence( x ,Q 2 ) = dy-Gu (y)U valenee (~,Q 2 ) (7) 

Jx y u y 

_ _i_ f^~X — X 

xd valence(. x >Q 2 ) = dy-G D_ (y)dvalenee (-, Q 2 ) (8) 
Jx y 7T+ d y 

xdZalenceix^Q 2 ) = I d y-G^_ (y)dvalenee (-,Q 2 ) (9) 

Jx y t D y 

XUvalence( x ,Q 2 ) = [ dy-G Q_ (y)U valenee (~,Q 2 ) (10) 

Jx y 7T- u y 

As for the sea quark distribution , we will take the example of ir + . ir + has 
two valons, U and D and each contributes to the sea and gluon content of 
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Figure 1: Comparison of the pion valence distribution, u^, alence (x), result from the valon 
model calculation (solid line), SMRS, Ref. [13] (dashed line), GRS, Ref. [14] (dotted line), 
and the data from E615 [2] at Q 2 = 25GeV 2 . 

pion as: 

Q 2 ) = f'dy-Gu (y) qw (-,Q 2 )+ C dy-Go(y)q^(-,Q 2 ) 
j x y n+ y jx y *+ y 

(11) 

Evaluation of these convoluted integrals requires us to determine G vai on (y) 

7T 

or, alternatively, finding /i m and v, m . Since the two valons in the pion, 
apart from flavor, cannot be distinguished and since the masses of the U 
and D valons are the same, therefore their average momentum also must 
be the same. This can be achieved only if [i m = v m , leaving us with only 
one parameter. We use the valence distribution data in ir~ at Q 2 = 25 
GeV 2 from the Drell-Yan experiment of E615 collaboration [2] to find this 
parameter. The fit to the data of Ref. [2] is shown in Fig. (1), which is 
obtained by taking the starting scale, Qq = 0.47 GeV 2 with Aqcd = 0.22 
GeV. 

The goodness of the fit is checked by x 2 minimization procedure. We 
find that the x 2 P er degree of freedom is 1.2 and that fi m = v m = 0.01. As a 
comparison, in Fig.(l) we have also shown the results from the determination 
of [13] as dashed line and those of [14] as dotted line, both at Q 2 = 25 
GeV 2 . In fitting the data we have also allowed both Aqcd and Qq to be 
free parameters and obtained an equally good description of the data with 
Ql = 0.35 GeV 2 and Aqcd = 0.175 GeV resulting in fi m = v m = 0.03. 
However, it seems that Aqcd = 0.175 is too low for four active flavor. As 
such, we choose the former against the latter. With the Determination of 
these parameters, it is important to check the valence quark sum rule in 
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pion. Using the explicit form of Eq. (2), i.e. zq™i°n Ce (z = x/y,Q 2 ), from 
the appendix of Ref. [11], we find that each of the integrals in Eqs. (7- 
10) gives 0.9994, 1.002, and 1.008 at Q 2 = 3,10,20 GeV 2 , respectively; an 
excellent confirmation of the valence quark sum rule. The first two moments 
of pion valence quark distribution is also calculated at Q 2 = 49 GeV 2 for 
the purpose of comparison with the findings of Ref. [13]. We find that 
2 <C XQ valence >= 0.378 and 2 < X Q valence >= 0.151. This is to be compared 
with 0.40"± 02 and 0.16 ± 0.01 of Ref." [13], respectively. 

From Eq. (5) we can infer some knowledge about the charge and matter 
distributions of the pion [8] . The longitudinal momentum space is related 
to the coordinate space by a Fourier Transform and a boost to infinite- 
momentum frame. Since the valon structure originates from QCD virtual 
processes, which are flavor independent, the matter and the charge densities 
in the valon ought to be flavor independent. Hence, the charge density in, 
say, 7T+ is 

P q {y) = \Gu{y) + \Gj, (y), (12) 

whereas for the matter density distribution we assume that it is proportional 
to the total valon distribution, i.e. , 

p m (y) = \{Gu(y) + G^(y)). (13) 

The integral of both quantities in Eqs. (12) and (13) are equal to one. 
Let us write the explicit form of G vaion : 

"in - 

Gu_ = G r> = 1.020(1 - y) om y om . (14) 

It is evident that this function is very broad in momentum space. This 
feature is expected, for, it indicates that the valons are tightly bound. This is 
also a reflection that the pion is much lighter than the mass of its constituent 
quarks. The parameters, /x m and v m obtained here are slightly different from 
those used in Ref. [4,11,12] and are significantly different from those quoted 
in [15]. In [4,11] the values /i m = 0.01 and v m = 0.06 were used, while the 
determination of Ref. [12] is [i m = u m = 0, thus, there is no significant 
differences among this work and those of Ref. [4,11,12]. In Ref. [15] the 
pion cloud model in conjunction with the valon model is used to calculate 
the pion structure. They find that /j, m = 0.044 and v, m = 0.372. In this 
work and in Ref. [4] the parton distribution in a valon is derived from 
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QCD alone, with no phenomenological assumption. However, in Ref. [4], 
the focus was on the low x behavior of the pion structure function, whereas 
here we have used data at rather large x (x > 0.2) region to determine the 
valon distribution in the pion. In Fig. (2) we present the full pion structure 
function at Q 2 = 7 and Q 2 = 15 GeV 2 at small x, (x = [1CT 4 , 1CT 2 ]), region 
and compared the results with those obtained from Sutton, et al. [13] and 
Gluck, Reya, and Schienbein [14] parameterizations. In Fig. (2) two sets 
of data points are shown, which correspond to the two different methods of 
normalization used by ZEUS collaboration. While the results from Ref. [14] 
agrees with the additive quark model normalization, the parameterization 
of Ref. [13] is qualitatively closer to the effective one-pion-exchange model 
normalization. Note that although the parameterization of Ref. [14] provides 
a good description of the additive quark model normalization of the data, it 
fails to describe the large x data of Ref. [2] as is apparent from Fig. (1). As 
can be seen from the figure, the valon model results are in good agreement 
with the pion flux normalization of the data. It is also interesting to note 
that in the valon model a simple relationship holds rather well between 
proton and pion structure functions, namely, Fg = kF%, with k ~ 0.37 [4]. 
This observation has also been made by the ZEUS Collaboration [3]. The 
ZEUS group will soon release two new measurement: one is photoproduction 
study from which the pion trajectory can be determined; and the second is a 
measure of the exponential Pt slopes in deep inelastic scattering that can be 
used to limit the choice for the form factor F(t) [16]. These measurements 
should help to resolve the normalization issue. 



3.2 B. Kaon 

The treatment of the kaon structure function is similar to that of pion, 
except that we need to determine valon distribution in the kaon. We will 
concentrate on K~ , since there are some data [17] which provide information 
about the valence distribution in K~. In Refs. [8, 11, 12] it is stated that 
the general form of the exclusive valon distribution in a meson is as follows 

GWyi,y 2 ) = \P(jik + + i)]~V fc 2/2 fc %i + 2/2-1) (15) 

Integrating over either of yi will give the individual valon distribution in 
the meson, as in Eq. (5). So, we can write a similar equation for the 
valon distributions in K~ , except that pL m and v m are different and we label 
them as and v^, respectively. One way to determine these parameters 
is fitting G vaion (y) to some experimental data. Unfortunately, there is 
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Figure 2: Pion structure function at Q 2 = 7GeV 2 (Left) and at Q 2 = 15GeV 2 (Right). 
The diamonds and squares are pion flux and additive quark model normalization of the 
data [3], respectively. The solid line represents the calculated results from the valon model. 
The bdashed line is the result from SMRS [13] and the dotted line corresponds to GRS 
[14] determination. 



no experimental data on the valon distribution of any hadron, and from 
the theoretical point of view, the function G vaion (y) cannot be evaluated 

hadron 

accurately. However, there are some limited data on the ratio 

R=^- (16) 

xu w - 

at large x values [17], and can be used to fit equations of the form (7-10) 
to get Hk and v k . Maintaining the same values for Qq and A as in the 
case of pion, we have fit the data of Ref. [17] and obtained n\. = 0.13 and 
Vk = 0.28 with a x 2 = 0.646 per degrees of freedom. The fit is shown in 
Fig. (3). In our fit to this data set, both valence and sea u are included, 
although the data points are at rather large x and hence, the sea quark 
contribution to R is marginal. In the alternative, we can be guided by 
making a simple phenomenological assumption as follows [12]. Since K~ 
consists of two valons with different masses. It is obvious from Eq. (15) 
that the average momentum fractions of the light valon, yi an d the heavy 
valon, y 2 , are yi = {nk+t) / {Hk+Vk) andy 2 = (ffc+l)/(^+^fc), respectively. 
Thus, letting the ratio of the momenta be equal to ratios of their masses, 
we get 

yi = m u „ 300 = ^ + 1 = Q 6 {m 

y 2 m s 500 u k + 1 
In this way v k is restricted by 

v k = (ji k + 0.4)/0.6, (18) 
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Figure 3: The ratio R as a function of x as measured by [17]. The curve is obtained 
from the model by fitting the data at Q 2 = 25 GeV 2 . 
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Figure 4: The ratio R as a function of x as measured by [17]. the curve is obtained with 
fit = -0.35 at Q 2 = 25 GeV 2 . 

hence, leaving us with only one parameter, fi^, to be determined. A one 
parameter fit to the data of Ref. [17] is highly unsatisfactory, with fx^ = 
—0.35 as can be seen in Fig. (4). Authors of Ref. [12] have used inclusive 
distribution of K + — > ir + and K + — > ir~ and found jjl^ = 1 and ^ = 2 
which also do not support the data of Ref. [17]. An earlier QCD calculation 
[18] of the ratio, which takes into account the difference in quark mass, 
however agrees with the data. Nevertheless, it seems that R being around 
0.5 at x = 0.95 is too high. Considering the large error bars in the data, 
the accuracy of the data, at least, at large x is suspect. One possible way 
of determining the valence quark density in kaon would be the difference of 
two cross sections in K + p — > / Lt + / u _ + X and K~p — > p~ + X processes. 

While the determination of the valon distribution in kaon remains un- 
certain, we will proceed with the values obtained for and \ik without the 



9 



0.2 0.4 0.6 0.8 



Figure 5: The valence quark distribution, xu(x), in n (solid curve) and in K (dashed 
curve) . The curves are obtained from the model at Q 2 = 25 GeV 2 . 

restriction of Eq. (18), namely, ^ = 0.13 and = 0.28 and arrive at the 
following valon distributions in K~; 

G k D (y) = 1.4768y°- 13 (l - yf 2S , Gf(y) = 1.4768y a28 (l - y) 013 

(19) 

With these relations the average momentum fraction carried by U and S 
valons are 0.469 and 0.531, respectively. Using these valon distributions 
along with Eqs. (7-10) for K~, the valence quark distributions for tt~ and 
k~ are calculated at Q 2 = 25GeV 2 and shown in Fig. (5). 

It should be noted that in Drell-Yan processes, obtaining experimental 
information on the strange valence and sea quark distribution in kaon is 
not practical. Because those components in kaon only contribute to the 
total cross section of the Drell-Yan processes through valence-sea and s — s 
annihilation and these contributions being small makes it difficult to separate 
them. The valon model, on the other hand, provides valuable information 
about these components. Figures (6) shows the strange valence and sea 
quark distributions in K~. Both are calculated at Q 2 = 25 GeV 2 from the 
model. For the purpose of comparison, in Figure (7) we show sea quark 
distributions in a valon. The sea quark distribution is calculated originally 
for u and d type partons. The strange sea quark distributions are generally 
smaller than the light sea quark distributions so that in proton 2s/ (u + d) ~ 
0.5[19] Nevertheless, by the time x ~ 10 -4 s w u. 
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Figure 6: Strange quark distribution in K~ at . Q 2 — 25 GeV 2 . Left figure is the 
valence component and the right figure is that of strange sea distribution. 
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Figure 7: Sea parton distributions in a valon at Q 2 = 25 GeV 2 . 
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4 Conclusion 



We have used the notion of the valon model to determine the structure of 
mesons. Since the partonic structure of a valon is already known, the struc- 
ture function of any meson needs only two parameters to be completely 
determined. The valon model provides information about the sea quark dis- 
tribution of mesons that are out of reach experimentally. While the pion 
structure is determined fairly accurately, the structure of kaon remains un- 
certain, due to the lack of accurate data or incompatibility of data sets from 
different experiments. 
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